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Abstract 



I consider gravity induced over a smooth (finite thickness) sohton. Gravi- 
ton kinetic term is coupled to bulk scalar that develops solitonic vacuum 
expectation value. Couplings of Kaluza-Klein modes to soliton-localized mat- 
ter are suppressed, giving rise to crossover distance Tc = Mp/M^ between 
4D and 5D behavior. This system can be viewed as a finite thickness brane 
regularization of the model of Dvali, Gabadadze and Porrati. 
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I. INTRODUCTION 



Brane worlds theories with large or infinite extra dimensions recently provided insight 
to a number of problems of high energy physics and possible relations between them. Those 
models can emerge as a part of the fundamental higher dimensional theory described in 
terms of strings. Even if extra dimensions are not realized in nature, they are valuable 
testing ground for new ideas in particle physics and gravity. 

For some time, it was believed that the extra dimensions can exist only if they are 
compact or have a finite volume. The problem with the infinite volume extra dimensions 
was how to obtain a four dimensional force of gravity between objects located at the brane. 
In work [1], it was shown that the four dimensional Newton's force arises on the brane due 
to an induced kinetic term. The model of Dvali, Gabadadze and Porrati [1] is described by 
the action 

S^Ml j d'^xdy^ln^ + Mj>J d'^xdy5{y)^\Ri. (1.1) 

Here y is the coordinate of an extra dimension, TZ^ is 5D curvature, Gmn 5D metric, 
is 4D curvature constructed from metric induced on the brane g^j^i, = GuNiy — 0)S^S^ , 
M* is the fundamental scale of gravity and Mp scale that characterizes brane rigidity. In 
order to obtain the correct vahic for Newton's constant on the brane, Mp has to be taken 
~ lO^^GcV. The brane is represented with the delta function i.e. thickness of the brane is 
taken to be zero. A related model with one infinite extra dimensions was presented and 
studied in [2]. Phenomenology of the action (1.1) was studied in [3]. There it was shown 
that the scale of the bulk gravity can be as low as M* ~ 10~^eV. Further properties of the 
model, cosmological evolution, new manifestations at astronomical scales and string theory 
realizations were examined in [4-6]. 

In these notes, I study a similar system of induced (enhanced) gravity over a smooth 
soliton 

S^Ml j d^xdy^\n,{l + ^x\Mpy)). (1.2) 

Here x(Mp7/) is a localized function representing a smooth tension-less soliton. I find that 
at the hnearized level, features of (1.1) and (1.2) are the same. In short, four dimensional 
gravity at distance r < Tc = Mp/Ml is mediated by a metastable massless resonance. Width 
of the resonance is ~ r~^. At distances r > gravity becomes five dimensional. Therefore 
I conclude that (1.2) is a finite thickness brane regularization of (1.1). 

Through the paper I will work with the simpler 'scalar gravity' model that captures main 
features of gravity. A Variation of the strength of the kinetic term can be a consequence of 
dilaton-like coupling to bulk scalara (x) that exhibits sohtonic vev. I neglect the gravitational 
effect of the soliton, which is equivalent to fine-tunning its tension to zero. In the next section 
I study properties of the spectrum of Lagrangian in which the strength of kinetic term is a 
function of the extra coordinate. Results arc applicable for scalar and gauge bulk field as 
well as for linearized gravity. In section 3, I investigate the spectrum of Kaluza Klein modes 
in system (1.2) and their coupling to soliton- localized matter. I find the crossover distance 
= Ml /Ml, both from the suppression of couplings of KK modes to matter localized on 
soliton and from an equivalent problem of tunneling through the potential barrier. I find 
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no new scale ~ \Jfc/o. (where 1/a is the brane thickness), contrary to results of [7], where 
the different regularization scheme was considered. I conclude that the system (1.2) can be 
viewed as a finite thickness brane regularization of action (1.1). In the appendix 1 outline 
the procedure for finding the spectrum of fermionic modes when the kinetic term of a bulk 
fermion field is varying through the extra dimension. 



II. LOCALIZATION IN MODEL WITH NON-HOMOGENEOUS KINETIC TERM 

In this section I investigate general properties of a model in which the strength of a 
kinetic term for a bulk field is not homogeneous through the infinite extra dimension. Let 
us take that the bulk field $(a;^, y) is coupled through the kinetic term to a field x(a;^, y) 

C = f{x)dA^d^^ + dAXd\ - V{x). (2.1) 

Function f{y) characterizes the coupling of two fields. Field x(a:^,|/) has potential V{x) ^"^^ 
can obtain an expectation value in the form of a soliton x — ^ Xciiy)- Then, the Lagrangian 
for the field $ has the following form 

L = f{y)dA^d''^. (2.2) 

Field $ can represent bulk scalar, gauge field or graviton in linearized approximation. 

To find the four dimensional spectrum of (2.2) we decompose field $ into Kaluza Klein 
modes (^{x^,y) = $m(2/)cm(^/i)) where — —rn^am- The differential equation for 
wave functions $m(2/) is 

fiy)dl<^m + dyf{y)dy^m + fiyW^m = 0. (2.3) 

Taking equations (2.3) for two different eigenvalues m, n, cross-multiplying them with eigen- 
functions, subtracting and integrating over the length of extra dimension we get the orthog- 
onality relation 



(2.4) 

where am is the normalization constant. Combining (2.3) and (2.4) we get 

Ifiy)^n^mdy ^'-'^ 

For non- negative f{y), both the denominator and numerator of this expression are positive, 
which means that there are no ghosts, tachions or ghost tachions in effective four dimensional 
action obtained from (2.2). We can introduce additional coupling of field $ to a 'charge' 
density '^{y)p{x) localized at the soliton (centered at y = 0). In 4D effective theory the 
coefficient of a coupling is determined by the convolution of the matter profile ^(y) and the 
wave function (an average of wave function over the soliton- localized matter) 

$„(0) = J ^{y)^Mdy. (2.6) 
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Effective action for modes then reads 



^-\Y. (d.^md^'^m - m'al) + -^^^a^p{xl (2.7) 



where is the fundamental scale. Zero mode ctq will be normalizable for choices of f{y) 
whose integral is finite, while it will be a single mode in continuum for functions f{y) that 
are nonzero at y — > ±00. This feature distinguishes between the cases of localization 
and quasi-localization. In order to find the spectrum and determine localization properties 
of KK modes, one has to specify the function f{y). To better understand the connection 
between the shape of f{y) and localization properties, it is useful to transform (2.3) to 
related Schroedinger problem 

^m = ^, -dl<Pm + V{y)(l)m^m^<f>m, (2.8) 

where the localizing potential V has the form 

Whether this potential will behave as an attractive well or a repulsive barrier depends on 
the balance between the magnitude of the second and first derivative of the profile f{y). 

A. f{y) = e-y''^' 

Now we consider two choices of function f(y) that lead to the localization of part or all 
of the spectrum of KK excitations. First example is a Gaussian profile f{y) = " . Here, 
a is the scale that determines the brane width. This profile will give rise to a localized tower 
of KK modes with masses m ~ y/n, where n is the integer. This was used as a model of 
localization of gauge fields in [8]. Here we repeat it because we will use it as a starting point 
for the investigation of soliton induced gravity in the next section. 

Potential (2.9) that results from this profile is the potential of a simple harmonic oscillator 
shifted by a constant so that the zero point energy vanishes 

y(y) = a2(aV-l)- (2-10) 
Normalized wave functions (2.8) are 

(f)Jy) = , e-^y^'^'/^HJya), (2.11) 

where the spectrum of masses goes like 

mn^aV2n, n = 0,1,2,... (2.12) 

The interaction of those modes with some charge located at the brane is determined by the 
value of the wave function at the position of the charge. For example, coupling of a charge 
at y = is 



4 



I J van. 

for n even, while it vanishes for odd n. The vanishing of couphng for odd values of n is also 
true for a smooth brane that locahze charge p. The reason is that wave functions of modes 
with odd n are odd, while the matter localization profile ^'(y) is even and coupling (2.6) 
vanishes. 

B. fiy) = e-\y-\ 

Other example is the exponentially decaying profile f(y) = e~'^"'. This is presented as 
an example of a localized profile f{y) with a finite integral, that doesn't localize all of the 
KK modes (as one could naively expect). The Schroedinger potential for this profile is 

V{y)^^a'-a5{y). (2.14) 
Spectrum of modes consist of a single localized zero mode 

00 - V^e-«l^l/^ (2.15) 

and a continuum of scattering states starting at m > a/2. The potential between two charges 
qi and q2 located at y = due to exchange of field $ would change from a four dimensional 
law to a five dimensional law at distance r ~ 1/a. 

III. SMOOTH BRANE REGULARIZATION OF DGP MODEL 

In this section I establish the relation between the model of Dvali, Gabadadze and Porrati 
and action (1.2) at linearized level. Let us consider the five dimensional model with two 
scales. Fundamental scale determines the action of the bulk filed $. The potential 
for the other bulk field x is governed by the much higher brane scale a. Therefore, scale a 
determines the width of the brane, which is soliton of field x- Field x is coupled quadratically 
to a kinetic term of and the coupling is suppressed in powers of fundamental scale 

dA^d^^ + dAxd\-V{x)- (3.1) 

Let us take that the field x gets expectation value xci — ot^^^ exp(— (ya)^/2) so that the 
profile function (2.2) is 

Since we take that the sohton scale a is much bigger than the bulk scale M*, coefficient A 
is a huge number. The Schroedinger potential (2.9) for this profile is 
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V{y) 



-Aa 



,A{l-{yaY) + e^y^y{l-2{yaf) 



(3.3) 



Potential (3.3) is shown in Fig. 1 together with the potential of harmonic oscillator (2.10). 
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Figure 1: Potential (3.3) (solid line) together with the harmonic oscillator potential (2.10) (dashed 
line), y is in units of inverse a, the potential is in units of a?, while A = 10000. 

For small values of y, the potential coincides with the potential of a harmonic oscillator. 
After certain a distance ~ 1/a, the square well is 'cut' and the potential rapidly drops to 
zero. From this behavior we can qualitatively predict the spectrum of states. Lowest bound 
states of the harmonic oscillator will obtain width and will tunnel through the barrier. Since 
we are interested in the case of gravity we want to find the decay width of a zero mode that 
will mediate the four dimensional Newton's force on the brane. Higher modes have mass 
~ a and are irrelevant for low energy physics. Decay width in the WKB approximation for 
the zero mode is 



Fq ~ aexp —2 



(3.4) 



where t/i,t/2 are classical turning points and 1/a is the width of the well. For large values 
of A condition for WKB validity |V/(2\/)^/^| ^ 1 is well satisfied everywhere (except at 
classical turning points). The integral in (3.4) can be evaluated numerically for the large 
range of values A. The integral is independent of a and behaves as ~ ^o-^ss (power has 
negligible numerical error). Therefore, in the WKB approximation zero mode has a decay 
width 

(3.5) 



^0.976 ■ 



Wave functions of the continuous Kaluza-Klein spectrum can be found by solving the 
Schroedinger equation. We find that the couplings of modes (2.6) to matter localized on the 
soliton are suppressed. 

Here, we can explore the question of the universality of the couphng of filed $ (scalar 
gravity) to brane matter. If different types of matter localized on the soliton have different 
profiles ^(y), coupling (2.6) of gravity would be non- universal. This is sometimes argued 



6 



to be a problem ^ for theories with gravity induced over a smooth brane [9] . We adopt the 
following framework in which the violation of universality is suppressed by soliton thickness 
scale a. The potential for a scalar field that forms solitonic brane is governed by high mass 
scale a. Here we take that a ~ Mp. Mass of Standard Model fields is much lower than 
the cutoff a. We take that SM particles come from the zero mode of another bulk field, 
while the mass splitting is due to some higher order effect. Zero modes have profiles roughly 
of a soliton itself \E'o ~ aexp(— (ya)^). Potential well that localizes zero mode can have 
bound states whose masses would be ~ a. Indeed wave functions of those states ^'s.s. are 
different and lead to the different couphng (2.6). However, because of their large mass they 
are decoupled from low energy physics. Now let us return to the spectrum of potential (3.3). 
The Schrocdinger equation for (3.3) can be solved numerically. For simplicity, we take that 
the matter profile is \E'o ~ aexp(— (ya)^) and evaluate coupling (2.6) for the continuum 
of KK modes. The qualitative behavior of wave functions is as follows. Outside the well, 
wave functions are plane waves of frequency m. Inside the well, wave functions are roughly 
Gaussians with mass-dependent height and constant width (~ 1/a). Approximately, one 
can take for the coupling ^m(O) ~ y^\E'm(0). Actual coupling was evaluated for each mode 
numerically. Dependence of coupling (2.6) on the mass of KK mode is shown in Fig. 2 

. 1 I I I 
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Figure 2: Suppression of KK coupling for A = 10000. For m = coupling has value 1 (outside the 
graph range). Peaks are at positions of resonant states of harmonic oscillator m = 0, 2, \/8... . Log-log plot 
of the same graph (multiplied with VA) is shown in the inset. The solid line has coefficient -1 to show the 
suppression of coupling ~ 1/m. 

There is a peak located at the mass zero. This is the zero mass graviton resonance that 



^non-universality of coupling can in fact be a desirable property of models with infinite extra 
dimensions. It opens the possibility for solving the cosmological constant problem without spoiling 
general covariance of the resulting four dimensional theory. This was pointed out to me by Gia 
Dvali. 
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is responsible for the 4D force at short distances. This resonance starts at couphng 1 and 
decays as ~ 1/m (see the log-log plot). Dependence of the half-width of the peak at m = 
on the magnitude of A can be evaluated by numerical analysis. By sampling data for the 
large range of values A {10^ < A < 10^) we find 

~ ^0.997±0.005 - (3-6) 

Besides the metastable 'massless' graviton there is no other resonances until we get to 
higher states of the harmonic oscillator that have mass ~ a. In particular, nothing unusual 
is happening at scale ~ a/y/A contrary to the results of [7] where a different regularization 
was used (there advocated mass scale a/y/A would be on position -4.605 on the log-log graph 
in Fig.2). 

Newtonian potential due to an exchange of field $ between two masses at short distances 
is given by 



V(r) ^ [ 
Jo 



^^^-^ dm ^ — r < A a. 3.7 

M3 r M^Ar' ' ^ ' 



Here the integral is supported over the resonance of width ajA and height 1. At large 
distances r ^ Aja only the lightest modes (m < ajA) contribute to the exchange. Their 
coupling is unsuppressed ($^(0) ~ 1) and one can approximate 

Vir)^ ' 7)7' dm^—-, r^A/a. (3.8) 

This is the same behavior of Newtonian potential as in the model of Dvali, Gabadadze and 
Porrati. A weak four dimensional force is mimicked by the exchange of a continuum of 
massive KK gravitons at distances smaller than the crossover distance r < = A/ a. At 
larger distances, the force becomes five dimensional. Crossover distance in (3.1) coincides 
with the crossover distance in [1]. In order to reproduce the correct Newton's constant one 
has to take the scale a to be of the order of Planck scale Mp 

_A _ a^ _Ml 

At the linearized level, an exchange of massive spin two states gives 'wrong' tensor struc- 
ture of the propagator [10]. It was shown [5,11] that non-linear effects cure this problem in 
model [1]. I conjecture that a similar mechanism occurs in the model at hand. Investigation 
of that effect as well as the spectrum of induced gravity in codimensions two and three is 
currently in progress [12]. 
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APPENDIX A: METHOD FOR FERMIONS 



Here we outline the procedure for finding the KK spectrum for bulk fermion field with 
kinetic term that is not homogeneous in extra dimension. The action we consider is 

S^J d'Xf{y)m^Ad^^. (1.1) 

Dirac algebra in 5D {7"^, 7^} = 2g^^ can be satisfied with a choice of 4D gamma matrices 
7^ and 7(^5^) = -i7f4C) = 'h^- 

The equation of motion is independent of f{y) and is solved with the ansatz 

'^{p,y) = Ei9m{y)^T + hm{y)^^) 

m 

7Xer = K]?, l,p'^R=m^T, (1.2) 

where spinors and together form a 4D Dirac spinor -0^, i.e. |(1 — 7^)-^^ = |(1 + 

l^)4'm — Cr- Wave functions are 

gm{y) = Am sm{my) + Bm cos(my) 

hm{y) ^ AmCoa{my) - Bmsin{my). (1.3) 

Although wave functions arc the same as in the case /(y) = 1, four dimensional effective 
action will be different, since one integrates over the profile f{y)- Plugging (1.3) into (1.1) 
and integrating over y we obtain the following 4D action 

m,n \ ^ / 

where is an infinite dimensional matrix given by a Fourier transform of the profile 

fmn = j dyf{y) cos(m - n)y. (1.5) 

Clearly, finding the spectrum of KK fermions (as well as their couplings) is now reduced to 
a problem of diagonalization of infinite hcrmitian matrices K = fmn and M = {m-\-n)fmn/'2- 
Diagonalization can be done in steps, and we can write schematically 

t/jUpKi/j + i/j^Mi/j — > i/^jXiipi/ji + i/j^UMU^i/j — > 
ipHpiP + ^^t^^^^yt^. ^ ^(V;^^^^„ + miPm^m), (1.6) 

where we rotated fermion sequentially as ip — > Uip^ipi — \/\iip,ip — Vip. Orthogonal 
matrices [/, V are chosen to diagonalize the kinetic and mass term. The simplest example 
is f{y) = 1 in which case fmn = 5{m — n) and we recover ordinary KK tower of massive 
fermions. Another simple example is f{y) — 5{y) in that case fmn has 1 at every entry. 
All eigenvalues of that matrix are zero except one that is 2N + 1, where A'" is (infinite) 
number of positive modes in -0 (modes go from —m to m). Thus, only one mode will 
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have nonzero kinetic term. Matrix U that diagonahzes fmn in this case has entries 1 on 
the diagonal that goes from the lower left to upper right corner, entries -1 in first column 
and entries 1 in last row. while all other entries are zero. Matrix UMW has zeroes on 
diagonal. Thus, a single propagating mode has the mass zero. It would be interesting to 
find the spectrum of induced 4D fermions by diagonalization of the infinite matrix = 
/ dy{\ + 74exp(— (ya)^) cos(m — n)y and compare it to the results obtained by inducing 4D 
kinetic term on the brane [13,14]. 
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